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Abstract 

The Bragg, Leibfried, and modified Leibfried numbers are denned in the context 
of a theory of dislocation-mediated melting, and their values are determined from the 
properties of the dislocation ensemble at the melting temperature. The approximate 
numerical coincidence of the Bragg and modified Leibfried numbers is explained. 
The parameter K in the definition of the modified Leibfried number is shown to 
be the natural logarithm of the effective coordination number. Our analysis reveals 
that the Bragg number can be considered an elemental constant, in contrast to the 
Leibfried and modified Leibfried numbers. 
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In Leibfried's study of melting |1[ he observed that for a number of metals (specifically 
Al, Ag, Au, Cu, Pd and Pb) 

L = RT m /GV ~ 0.042, (1) 

where R is the gas constant, T m is the melting temperature, G is the ambient shear 
modulus, and V is the molar volume. The factor RT m is an approximation to the heat of 
fusion, L m , which follows from Richard's rule |§ 

L m /T m = AS m ~ 2.0 e.u. « R (2) 
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(1 e.u. = 1 cal/mol °K), since the numerical value of R is 1.987 e.u. Bragg ||, on the 
other hand, noted that 

B = L m /GV ~ 0.034 (3) 

for a few metals (specifically Al, Ag, Au, Co, Cu, Fe, Ni and Pb). 

Although L, the Leibfried number, and B, the Bragg number, seem to be intrinsically 
related, through Richard's rule, there is a relatively large discrepancy between their values. 
The importance of both numbers in the calculation of the size factors caused Gschneidner 
to reanalyze them in his review ||. First, he noticed that Richard's rule may not be 
a valid approximation. In fact, Stull and Sinke j|] had earlier used two different values 
for estimating heats of fusion, 2.3 e.u. for elements which crystallize as face-centered 
cubic (fee) or hexagonal close-packed (hep) structures, and 1.9 e.u. for elements which 
have body-centered cubic (bec) structure. Second, Gschneidner introduced the modified 
Leibfried number, L', which differs from the Leibfried number L in that the term RT m 
is replaced by the term KT m . The value of K depends on the crystal structure of the 
element just below its melting point: 2.29 e.u. for fee or hep metals, 1.76 e.u. for bec 
metals, etc. Finally, he concluded that the modified Leibfried and Bragg numbers do 
agree within uncertainties. 

In this Technical Note we offer an explanation for the numerical values of both V 
and B in a theory of dislocation- mediated melting. We obtain the numerical values of B, 
L and V from the properties of the dislocation ensemble at the melting point. We also 
explain the approximate numerical coincidence of the values of V and B. 

In our previous study of melting as a dislocation-mediated phase transition on a lattice 
|5|, |(|, we obtained two relations: 

\kGv ws . ( a 2 \ 

L m = jb 2 p(T m )RT m In z', (5) 

where p{T m ) is the critical dislocation density at melt, vws is the Wigner-Seitz volume, 
1/k = (1 — f/2) ± u/2 « 5/6, v being the Poisson ratio, b is the shortest perfect- 
dislocation Burgers vector, A = b 3 /vws ~ 4/3, and a = 2.9 accounts for nonlinear effects 
in a dislocation core. Also, z', the effective coordination number for a dislocation as a 
random walk on the lattice, satisfies z' < z, where z is the coordination number of the 
lattice. These relations hold to ~ 20% accuracy for those elements of the Periodic Table 
(more than half) for which sufficient data exist. 
It then follows that the Bragg number is 

We have found || that the critical dislocation density at melt is an elemental constant 
with an approximate numerical value of 2/36 2 for three-quarters of the Periodic Table. It 
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corresponds to the situation where, on average, half of the atoms are within a dislocation 
core at melt. Hence, the Bragg number must be an elemental constant, and its approxi- 
mate numerical value is, in view of Eq. (6) with the numerical values of the parameters 
quoted above, 

92 

B « _ « 0.037, (7) 

8n 

in agreement with Bragg's original estimate. We discuss this point in more detail below. 
Similarly, the Leibfried number is 

Xk , / a 2 \ 

ln TI^T^T • (8) 



87rlnz' \<ib 2 p(T n 

It depends, through biz', on the crystal structure of the solid phase of the element from 
which melting occurs, in agreement with Gschneidner's observation. It is, however, pos- 
sible to define modified Leibfried numbers that approximately coincide numerically with 
the Bragg number ||: 

where K may be defined in either of two ways: 

g = 26V(r)ln/ eu (w) 

A 



or 



K = lnz'e.u. (11) 



In the first case, 



k , n . . / a 2 \ e.u. 
£ ' = _^T.)ln( 3p ^j- (12) 

will approximately coincide numerically with B given in Eq. (6), since the numerical value 
of R in units of e.u. is very close to 2. In the second case, 

L'=-ln( f r (13) 
8vr \4b 2 p(T m )J R 1 ; 

will again approximately coincide numerically with B given in (6), since R ~ 2 e.u., 
b 2 p(T m ) « 2/3 and A sa 4/3. 

In the first case, the numerical values of K in e.u., according to Eq. (10), are 2.26 for 
fee and hep metals (for which A = \[2 and z = 12) and 1.99 for bec metals (for which 
A ~ 1.3 and z — 8), where we have taken b 2 p{T m ) = 2/3 and z' — z — 1 || ||. These values 
of K are consistent with those presented by Gschneidner, 2.29 and 1.76, respectively. In 
the second case, the corresponding values for z' — z — 1 are 2.40 for fee and hep metals 
and 1.95 for bec metals, and for z' = z — 2 these values are 2.30 for fee and hep metals and 
1.79 for bec metals. It is therefore seen that in the second case the choice z' = z — 2 leads 
to better agreement between the values of K calculated from Eq. (11) and those given 
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by Gschneidner. This is not sufficient, however, to conclude that the choice z' = z — 2 
is better than the choice z' = z — 1 used in || ||, since Gschneidner did not specify 
the criteria for his choice of the values of K. Equation (5) with z' = z — 1, for example, 
explains the Stull-Sinke rule, since it gives L m /T m = 2.25 e.u. for fee and hep metals, 
and 1.98 e.u. for bec metals, in good agreement with the values 2.3 e.u. and 1.9 e.u. used 
by Stull and Sinke. In any case, the uncertainty in the value of In z' associated with the 
choice between ln(z — 1) and \\i(z — 2) does not exceed 8%, which is well within the ~ 20% 
uncertainty of Eqs. (4), (5). 

In ref. ||, for three-quarters of the Periodic Table, we established that 

b 2 p(T m ) = 0.64 ±0.14. (14) 

It then follows from Eqs. (3)- (5) with z' = z— 1 and the numerical values of the parameters 
quoted above, that 

B = 0.0369 ± 0.0093, (15) 

where most of the error comes from uncertainty in the value of k. Hence, uncertainty in 
the numerical value of B is about 25%, which is sufficiently small that B can be considered 
an elemental constant. 

As follows from Eqs. (1), (4) with z' = z — 1, and Eq. (14), 

L = 0.0308 ± 0.0113, (16) 

where we have used z = 8 and 12 to fix the error bar. Also, from using this value of L, 
the modified Leibfried number defined by Eqs. (9) and (11) with z' = z — 1 is 

L' = 0.0337 ± 0.0119, (17) 

where, again, the error reflects the two values of z. These values of L and V are in 
agreement with the values obtained in ref. |2|] for almost all the elements of the Periodic 
Table: L = 0.0305 ± 0.0135, V = 0.0334 ± 0.0145. Although the values of B, L, and V 
do agree within errors, uncertainties associated with both L and V are about 35% which 
is so large that they cannot be considered elemental constants. 

To conclude, we have defined and evaluated the Bragg, Leibfried and modified Leibfried 
numbers in the theory of dislocation-mediated melting. We have explained the approx- 
imate numerical coincidence of the Bragg and modified Leibfried numbers. We have 
shown that the Bragg number can be considered an elemental constant, in contrast to the 
Leibfried and modified Leibfried numbers. In our analysis, the quantity In z' is identified 
with Gschneidner's constant K. The factor Inz' can only arise naturally in a theory of 
line-like defects. 
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